Superradiance in ultracold Rydberg gases 
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Experiments in dense, ultracold gases of rubidium Rydberg atoms show a considerable decrease of 
the radiative excited state lifetimes compared to dilute gases. This accelerated decay is explained by 
collective and cooperative effects, leading to superradiance. A novel formalism to calculate effective 
decay times in a dense Rydberg gas shows that for these atoms the decay into nearby levels increases 
by up to three orders of magnitude. Excellent agreement between theory and experiment follows 
from this treatment of Rydberg decay behavior. 



In recent years, ultracold atomic gases have been used 
to probe a variety of many-body phenomena such as 
Bose-Einstein condensation P, |3| and degenerate Fermi 
gases 0- In addition to collective effects due to parti- 
cle statistics, other manifestations of many-body physics 
have been explored, such as in slow-light experiments Q 
and in ultracold Rydberg gases (e.g. the diffusion of exci- 
tations through resonant collisions and the blockade 
mechanism jy]). Another important fundamental col- 
lective effect is superradiance, in virhich photon exchange 
between atoms modifies the behavior of the sample. In 
particular, cooperative effects due to virtual photon ex- 
change can lead to the formation of so called Dicke states 
0. These states are the symmetric superposition of all 
states with the same total excitation level for constant 
atom number N. Interest in Dicke states has grown re- 
cently because of their potential advantages in quantum 
information processing Q and their irnportance in the 
behavior of Bose-Einstein condensates 

In this Letter, we are interested in many-body physics 
involving photon exchange in an ultracold gas of Ry- 
dberg atoms. Because superradiance depends on the 
atomic density per cubic wavelength, and because radia- 
tive decay of Rydberg atoms takes place predominantly 
between the closely spaced upper levels, ultracold Ry- 
dberg gases are ideal systems to study superradiance. 
In fact, Rydberg atoms have many interesting proper- 
ties: their size can become comparable to the atomic 
separation, and they have huge dipole moments p ^ n^, 
where n is the principal quantum number of the Ryd- 
berg state. In addition, for long-wavelength transitions 
between neighboring states of high n the "cooperative 
parameter" C ~ AfX^/An'^ (where Af is the density of 
atoms, A is the transition wavelength), is large for Ryd- 
berg atoms, which means collective effects are much eas- 
ier to obtain than for ground-state atoms ^3 • This was 
confirmed in earlier experiments for Rydberg atoms at 
high [nl and low temperatures 0| . Note that these 
many-body effects may pose a limit on the measurement 
of lifetimes of Rydberg atoms and may cause un- 



desirable frequency shifts, for example in atomic clocks 

The source responsible for both virtual and real pho- 
ton exchange is the dipole-dipolc interaction. It gov- 
erns the build-up as well as the decay of coherence in 
a dense radiating sample. On the one hand, the virtual 
exchange of photons is responsible for the so-called ex- 
change interaction. Its strength is exemplified by the 
energy difference 2h.^l = p^/27reor''^ between the sym- 
metric and anti-symmetric single-excitation superposi- 
tion |±) = (|eg) ± \ge))/V^ of two atoms in their ground 
g or excited e states separated by r. On the other hand, 
the real photon exchange is responsible for dephasing of 
a dense gas and has the same dependence. The inter- 
play of both determines whether the decay speed-up in a 
dense inverted gas of two-level atoms is mostly incoherent 
(intensity proportional to atom number N , called "am- 
plified spontaneous emission", ASE) or coherent (oc N'^, 
called "superradiance" or "superfluorescence" ) . Experi- 
mentally, this difference can be seen in whether there is 
an initial build-up in the decay intensity, due to the N'^ 
dependence, or not. 

The difficulty of calculating effects including atom- 
atom cooperation relates to the intractably large number 
of interconnected degrees of freedom, even if just a few 
particles arc involved. To explore these collective effects, 
many new ideas, such as the quantum jump approach, 
were developed to treat superradiance [l^ [l^, llal • R-^- 
cently, we successfully incorporated cooperative effects 
into a novel formalism for optically dense media. The 
result is a two-atom master equation for superradi- 
ance 0, I20I EH . We apply our model of cooperative 
radiation build-up to explain the results of an experi- 
ment measuring rapid decay of an ultracold Rb Rydberg 
gas. 

The model, as used in 0, , is based on perturba- 
tion theory carried to second order in the strength of the 
exchange interaction. Thus, we can eliminate all field 
and most atomic degrees of freedom which results in an 
effective two-atom nonlinear equation of motion of the 
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where p is a two-atom density operator, u^^ ^ is the lower- 
ing (raising) operator of the ith atom, 7 the spontaneous 
emission rate, and contains the second order dipole- 
dipole interaction between atoms i and j. (First order 
effects lead to local field effects which don't play a role 
here ^3)- order to obtain this result, Gaussian (and 
therefore, classical) light field statistics are assumed, in 
line with the second order approximation. In addition, a 
Markov approximation is made which is justified if the co- 
herence time of the light fields is shorter than the atomic 
evolution Atomic collisions and center-of-mass mo- 
tion are neglected. 

The Fy operators can be calculated from Tij5{t — t') cx 
{{Ei{t)Ej{t')), where Ei denotes the quantum field at 
the location of atom i, and the cumulant {{AB)) = 
(AB) — (A) (B). Fy's contain both the virtual and real 
photon exchange, and can be calculated for different sys- 
tems. They can be expressed only as highly nonlinear and 
implicit functions of the atomic variables p (Eqs. Q). For 
small enough probe diameters d retardation effects can 
be neglected. This approximation is justified in our case 
because the time it takes for light to propagate through 
the sample (~ 10""'^° s) is significantly shorter than any 
other time in the system, in particular, the atomic build- 
up time. Note that sample-sizes less than the cubic wave- 
length, as needed in the Dicke model are not neces- 
sary. Thus we can set F^^ = F and Tij^i = F and simplify 
Eq. CEJ: 
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The upper-level population is pee, the inversion product 
ITT- = (Pee ~ Pgg)^i ^ud the two-atom non-diagonal cou- 
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The sample size g = ird/X is measured relative to the 
wavelength of the light. 

In our initial experiment, we have studied the decay of 
high-n states using a simple detection scheme with only 
limited state-specificity. First, Rb atoms were trapped 
and cooled to 100 pK. Next, they were selectively ex- 
cited by a pulsed UV laser to the 40p state. After a 
delay time r, all atoms in states with principal quantum 
numbers n > 27 were Stark ionized. The remaining ex- 
perimental details are the same as in |2^ . As depicted in 
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FIG. 1: Measured and calculated decay of the number of 
atoms in Rydberg states with n > 27 following excitation 
to n = 43p. The initial density of Rydberg atoms in the ex- 
periment is 5 X 10* cm~"^ The dots are experimental points, 
the solid line theoretical simulation. The fitting parameter in 
this calculation was the number of atoms present at the start 
of the measurement, i.e., 1400 Rydberg atoms at 6 ps. 

Fig. n the number of ultracold atoms in Rydberg states 
with n > 27 decays fast, an estimated 100 times faster 
than expected in vacuum j^^ . We find that this speed-up 
can be explained by the presence of superradiance and, 
on some transitions, ASE [Hill [23. (Alternative ex- 
planations for the strength of the spee d-up would include 
so-called avalanche plasma formation |2J| , where a large 
fraction of the initial Rydberg atoms would be ionized. 
However, we rule this out because we measure only 190 
free ions after a delay of 35 ps.) 

In what follows, we will show that Eqs. ^ lead to ex- 
cellent agreement with the experiment (see Fig.QJ. The 
density in the calculation is chosen to be the same as in 
the experiment, 5 x 10® cm^^^. The sample in the experi- 
ment is cigar shaped, thus enabling good mode selection 
(as in all superradiance experiments to date). In the cal- 
culation we make the approximation of having, for each 
transition, only one mode, and then use, for calculational 
ease, a spherical geometry with the same sample volume 
as in the experiment. 

The simulations presented here assume Rb atoms in 
the initial state 40p. In Fig. |21 we show the decay from 
40p into ns. In Fig.|2| the effective decaytimes are com- 
pared for a dense gas and a vacuum (cf . |25| ) . In vacuum, 
,the effective decay time t^s is the inverse of the Einstein 
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FIG. 2: Calculated output intensity as function of time for a 
sample with density 5 x 10^cm~^ for the transition from state 
40p to 39s, 37s, and 6s, respectively. The initial increase 
in intensity over time is the sign for superradiance, i.e., the 
decay into 39s and 37s qualifies as superradiant, whereas the 
decay into 6s does not. (The curve for 6s is shown on a 100 
times faster timescale to show the decay.) 



A-coefiicient. Clearly, in a vacuum the transition into 
the states with lowest n is fastest, and therefore decay 
into these channels is by far the most likely. But this 
tendency is reversed dramatically in dense gases: the ef- 
fective decay time for each transition is shorter by up to 
three orders of magnitude than that in a vacuum or in 
dilute gases. Since the collective and cooperative effects 
responsible for this speed-up depend only on the density 
relative to the wavelength cubed, the acceleration of the 
decay is obviously stronger for longer wavelengths. Fig- 
ure 13 and the quantitative form of the increase in decay 
for higher densities, particularly for low frequencies, are 
one of the main results presented in this letter. 



must be proportional to the negative time derivative of 
the upper state population. (We neglect here all time- 
delay effects, resulting in an instantaneous intensity at 
time t ~ 0.) In this graph, our (somewhat arbitrary) 
distinction between ASE and superradiance can be seen: 
An initially positive slope of intensity over time, as seen 
for 40p — > 39s is associated with superradiance, whereas 
a monotonically decreasing intensity, as seen for 40p — > 
6s means ASE. It is important to emphasize here again 
that in reality there is no sharp boundary as there are 
coherent and incoherent elements mixed in all decays, 
thus making the transition between the two cases very 
smooth. 

In order to get a general overview of which combina- 
tion of parameters leads to superradiance, we created a 
map in the C-g parameter space with relative density 
or cooperative parameter C and relative size g. Figure 
0] shows the numerically determined border, as defined 
above, between superradiant and ASE behavior. The se- 
lective decay from the 40p Rydberg state of Rb into all 
possible lower ns states is added to the map. We see that 
superradiant behavior is expected for decay into levels 
with n> 22. 
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FIG. 3: Decay times from 40P to various nS states (o) in a 
dense gas (A/" = 5 * 10* cm"^) and in vacuum 

In Fig. 121 the intensity of some selective decays is shown 
over time. Because of energy conservation, the intensity 



FIG. 4: Map of critical parameters of C and g (o). Above 
the critical curve (shaded area) are the parameters for which 
superradiance happens. Also shown are the C and q for the 
decay to ns states from 40p state (-I-). Density of atoms is 
the same as above. 

We discuss now the calculation shown in Fig.^ Theo- 
retically, we can calculate, from the decay times as pre- 
sented in Fig.O the lifetime of 40p (and the lower states) 
directly, 

1 'ST- 1 
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and we find rtotai ~ 5fis. This is to be compared to a 
'''total ~ 210/J,s for dilute gas or vacuum. The experiment, 
however, cannot measure this time directly but only the 
total lifetime of all states with n > 27. In order to com- 
pare our theoretical method with the experiment we sim- 
ulate a cascade from 40p via all intermediate states down 
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to n < 27, using the decay times in Fig.|21 and analogous 
times for the p, d, /, etc. states with 5 < n < 40. This 
procedure is approximated by using only the two fastest 
channels out of each state. Numerically, we can compare 
this result with one that uses one channel more per state 
and find only small changes of 1-10%. The result is de- 
picted in the strong black curve in Fig.^ which shows 
excellent agreement with the experiment. 

Future experiments with improved state-selective de- 
tection will allow direct comparison to the single-lifetime 
calculations. 

In this article, we have discussed the possibility of 
superradiant decay in cold gases of Rydberg atoms at 
densities of 10^ - 10^ cm""^. Superradiance occurs be- 
cause lower-frequency decays are increasingly more likely 
to happen in denser gases, and they contribute most 
to cooperative behavior. Level shifts due to atomic in- 
teractions may inhibit superradiance at higher densities 
and/or higher n. This could explain why superradiance 
is not routinely seen. 

We have neglected black body radiation, since it is im- 
portant for superradiance only for the initiation of the 
radiation process and only if N ^ n^, where ub is the 
average number of black body photons per mode at the 
frequency of transitions |2(il |. In addition, the possibility 
of mode competition and interference between different 
decay channels is neglected for simplification. In future 
work, the effects of geometry, in particular the aspect 
ratio of the sample, should be taken into account. In 
practice, only elongated samples are used to show super- 
radiance [27I 1. 

Using our calculation we were able to obtain close 
agreement with observed signatures of superradiance in- 
cluding the effects of dissipation and the unique temporal 
build-up of a sharp flash of radiation. Moreover, our new 
formalism allows for easy incorporation of more compli- 
cated level structures, additional fields, and polarization 
effects. 

In summary, recent experiments measuring the decay 
of ultracold Rb Rydberg atoms find rates much faster 
than that of atoms in dilute gases. These results are 
consistent with superradiant behavior in the framework 
of our model. 
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